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Abstract 

Mg and Mg + concentration fields in the upper mesosphere/lower thermosphere (UMLT) 
region are retrieved from SCIAMACHY/Envisat limb measurements of Mg and Mg + 
dayglow emissions using a 2 -D tomographic retrieval approach. The time series of 
5 monthly means of Mg and Mg + for number density as well as vertical column density in 
different latitudinal regions are shown. Data from the limb mesosphere-thermosphere 
mode of SCIAMACHY/Envisat are used, which covers the 50 km to 150 km altitude re- 
gion with a vertical sampling of 3.3 km and a highest latitude of 82°. The high latitudes 
are not covered in the winter months, because there is no dayglow emission during 
io polar night. The measurements were performed every 14 days from mid-2008 until 
April 2012. Mg profiles show a peak at around 90 km altitude with a density between 
750 cm -3 and 2000 cm -3 . Mg does not show strong seasonal variation at mid-latitudes. 
The Mg + peak occurs 5-1 5 km above the neutral Mg peak at 95-1 05 km. Furthermore, 
the ions show a significant seasonal cycle with a summer maximum in both hemi- 
i 5 spheres at mid- and high-latitudes. The strongest seasonal variations of the ions are 
observed at mid-latitudes between 20-40° and densities at the peak altitude range from 
500cm -3 to 6000cm -3 . The peak altitude of the ions shows a latitudinal dependence 
with a maximum at mid-latitudes that is up to 1 0 km higher than the peak altitude at the 
equator. 

20 The SCIAMACHY measurements are compared to other measurements and 
WACCM model results. In contrast to the SCIAMACHY results, the WACCM results 
show a strong seasonal variability for Mg with a winter maximum, which is not observ- 
able by SCIAMACHY, and globally higher peak densities. Although the peak densities 
do not agree the vertical column densities agree, since SCIAMACHY results show 
25 a wider vertical profile. The agreement of SCIAMACHY and WACCM results is much 
better for Mg + , showing the same seasonality and similar peak densities. However, 
there are the following minor differences: there is no latitudinal dependence of the 
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peak altitude for WACOM and the density maximum, passing the equatorial region 
during equinox conditions, is not reduced as for SCIAMACHY. 


1 Introduction 

The amount of meteoric mass deposited on the earth upper atmosphere every day is 
estimated to range from 2 to 300 tons, (see, e.g. Table 1 of Plane, 2012, and references 
therein). 

Since the number of extraterrestrial particles encountering the Earth’s atmosphere 
is approximately inversely proportional, in a log-log scale, to the particle’s size, the bulk 
of the daily meteoric mass flux is contributed mainly by particles in the 1-100 microm- 
eter (in radius) size range (Ceplecha et al., 1998; Flynn, 2002; Mathews et al., 2001; 
Nesvorny et al., 2010, 2011a, b). Furthermore, most of the daily mass input originates 
from the sporadic meteor background, which is the interplanetary dust forming the Zo- 
diacal dust cloud and which orbits are so evolved that they cannot be traced back to 
their original parent body (Brown and Jones, 1995; Janches et al., 2006; Fentzke and 
Janches, 2008). The constant bombardment of particles from the sporadic background 
contributes far more input that meteor showers (e.g. Murad and Williams, 2002), the 
latter being a time limited (days to several weeks) enhancement of particles entering 
the atmosphere in parallel trajectories appearing to radiate from one point in the sky 
(Jenniskens, 2006). When meteoroids enter the atmosphere at geocentric speeds be- 
tween 11 to 72kms _1 they collide with air molecules, leading to frictional heating and 
deceleration. The meteoroids, consisting partly of metals, melt on their surfaces and 
metals are ablated. The ablation process depends on several factors, like the velocity 
and entry angle of the meteoroid, the boiling point of the different particles constituents 
and the size of the meteoroid (see, e.g. McNeil et al., 1 998 and Vondrak et al., 2008 for 
more details). For the micronsize particles of interest to this work, the ablation process 
occurs between 80 and 125 km. (Vondrak et al., 2008; Janches et al., 2009; Marsh 
et al., 2013a; Feng et al., 2013) resulting in the deposition of metallic atoms such as 
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Sodium, Iron, Potassium and Magnesium in the Mesosphere and Lower Thermosphere 
(MLT). The observation technique used in this study measures the density of the metal 
atom and ion layers present in the MLT that are formed with the ablated material from 
the meteoroids. Once ablated the meteoric metals undergo several chemical reactions 
5 (see e.g. Plane, 2003). A fraction of the ablated metals in the ablation region (80 to 
105 km) are already in the ionized form. For the case of Mg + this fraction is estimated 
to be negligibly small (Vondrak et al., 2008, Fig. 14.). Mg + is formed by charge ex- 
change with the main ion constituents at this altitude NO + and 0 2 . The dominant loss 
process for Mg is the reaction with 0 3 , which leads to stable oxides, hydroxides and 
io carbonates as reservoir species. Mg + is lost by reactions with 0 3 but also by reac- 
tions with N 2 and other trace gases (C0 2 , H 2 0) (see, e.g. Fig. 9 in Plane and Whalley, 
2012 for a schematic diagram of the most important reactions). The metal molecules 
condense and form meteoric smoke particles (e.g. Hunten et al., 1980; Kalashnikova 
et al., 2000; Saunders and Plane, 2006). These meteoric smoke particles can act as 
is condensation nuclei. Heterogeneous condensation of clouds plays an important role in 
the formation of polar mesospheric clouds in the polar summer mesopause region and 
also in the formation of polar stratospheric clouds (PSCs) in the polar winter strato- 
sphere. PSCs have an important role in the ozone chemistry and the formation of the 
ozone hole. In a dynamical equilibrium of meteoroid input and the above mentioned 
20 loss reactions, the neutral Mg layer is formed between 85 and 90 km altitude and the 
Mg + layer is formed slightly above the Mg layer (see, e.g. Plane and Helmer, 1995; 
McNeil et al., 1998 and Plane and Whalley, 2012). 

The MLT is not readily probed by in-situ instrumentation, as balloons fly too low, and 
satellites have to fly higher because of the atmospheric drag that reduces the satellite’s 
25 operation time. As a result, only rockets carrying in situ instrumentation are able to ac- 
cess this region. These are, however, expensive and rare. Fortunately metal atoms and 
ions have very strong absorption coefficients and oscillator strengths (see, e.g. Kramida 
et al., 2012). As a result, although their densities even in the peak region are only a few 
thousand particles per cubic centimeter, they are very strong emitters of resonance 
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fluorescence. Mg and Mg + spectral lines are observed at wavelengths below 290 nm. 
This wavelengths region is strongly affected by ozone absorption in the stratosphere, 
so that observations from ground are not possible. For this reason the knowledge of the 
Mg and Mg + content in the MLT is poor. In this study the dayglow emission of Mg and 
Mg + , detected by the satellite experiment SCIAMACHY (Scanning Imaging Absorption 
Spectrometer for Atmospheric CHartographY) on Envisat (Environmental Satellite) are 
used. The Mg and Mg + densities are retrieved by the mathematical inversion of these 
measurements using a radiative transfer model. Mg and Mg + have already been ob- 
tained from SCIAMACHY by Scharringhausen (2007) and based upon the results of 
this a new measurements mode for the purpose of better covering the peak region of 
Mg and Mg + has been established and used for this study. 

In the study these experimental data are compared to other Mg and Mg + measure- 
ments as well as simulations of the MLT using a combination of a Meteor Input Function 
model coupled to a meteor differential ablation model which provides the metal input in 
the MLT. The result of these models is then used as input in the Whole Atmosphere Cli- 
mate Community Model (WACCM) with a new chemistry model of Mg and Mg + (Plane 
and Whalley, 2012). WACCM is a general circulation model that incorporates interac- 
tive chemistry for both neutral and ion species. 

In Sect. 2 a short overview about the SCIAMACHY/Envisat satellite observations 
and the Mg and Mg + density retrieval is provided. In Sects. 3 and 4 the results for 
Mg and Mg + are presented and discussed. In Sect. 5 these findings are compared to 
other measurements. In Sect. 6 the WACCM-Mg model simulations of Mg and Mg + 
are introduced and the SCIAMACHY dataset is compared to the model results. Finally 
Sect. 7 summarizes the results obtained in this study. 


2 Instrument and algorithm 


SCIAMACHY, a satellite-borne grating spectrometer, was launched by ESA (European 
Space Agency) on Envisat on the 28 February 2002 (see, e.g. Burrows et al., 1995; 
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Bovensmann et al. , 1999). Envisat flies in a sun-synchronous orbit with a descending 
node (equatorial crossing time from north to south) at 10 : 00 a.m. The duration of one 
orbit is roughly 100 min (15 orbits per day). Note that the local time of SCIAMACHY 
limb observations changes with latitude. It is » 11 a.m. at 60° N and 9 a.m. at 60° S 
5 and, as shown in Fig. 1 , varies strongly in the vicinity of the poles. The highest latitude 
covered by SCIAMACHY limb observations of scattered light is 82°. 

Data from SCIAMACHY’s limb mesosphere and lower thermosphere (MLT) mode 
with tangent altitudes between 50 and 150 km in 30 consecutive 3.3 km steps are used 
in this study. Depending on latitude, consecutive limb measurements have a latitudi- 
io nal separation of up to 8 °, because nadir measurements are performed between two 
MLT-measurements. However, there is a small latitude offset between MLT measure- 
ments in two consecutive orbits. MLT-measurements were performed from mid 2008 
until the loss of contact to Envisat in April 2012. MLT-measurements were performed 
roughly every two weeks for 15 consecutive orbits, which corresponds to one day of 
15 consecutive measurements. Apart from a few days, where some orbits are missing 
in the data, all longitudes are equally covered (sampling every 24° of longitude). The 
two dimensional retrieval algorithm uses all measurements of one orbit and retrieves 
the metal densities from each emission line of the species on a latitude and altitude 
grid. For averaged data, a reference orbit for geolocations is used, which is formed 
20 from the dataset to average. Measurements with similar latitudes and local times as 
the measurements of the reference orbit are averaged. A comprehensive discussion of 
the algorithm and error estimation was made by Langowski et al. (2013). For Mg the 
spectral line at 285.2 nm is used. For Mg + the spectral lines at 279.6 nm and 280.4 nm 
are used. 

25 It has to be noted that the first measurements at the start of each orbit are con- 
taminated by solar stray light, which can be seen in the dark signal measurement at 
350 km tangent altitude. These measurements, as well as the next consecutive mea- 
surements, which do not show solar stray light in the dark signal, have been sorted 
out for this study. The strong solar stray light contaminated measurements actually 


V) 

O 

c 

C/) 

c n 

o' 

-o 

"0 

0) 

"O 

CD 


C/) 

o 

c 

(/) 

o' 

o 

"□ 

0) 

■O 

CD 


ACPD 

14, 1-49, 2014 


Global investigation 
of the Mg atom and 
ion layers 

M. Langowski et al. 


Title Page 


Abstract ■ Introduction 


Conclusions ■ References 


C/) 

o 

c 

V) 

c n 

o' 

o 

~0 

0) 

"O 

CD 


O) 

O 

c 

V) 

V) 

o' 

:d 

"0 

0) 

■O 

CD 


Tables 


Figures 



Full Screen / Esc 


Printer-friendly Version 


Interactive Discussion 


© CD 


7 

















show reasonable profiles. However, the measurement in between the stray light and 
not stray light contaminated measurements showed strongly different features than the 
neighboring measurements. At some days similar features for the remaining highest 
northern latitude measurements (after excluding the stray light measurements) still can 
5 be seen. This features are not dependent of the orbit phase anymore, but only appear 
at latitudes > 78° N. The actual grid for calculating the densities has a larger extend and 
the outermost latitude intervals have been cut off to reduce edge effects, which may, 
however, affect the neighboring latitudes. 


3 Seasonal variations of Mg 

io The seasonal variations in altitude and latitude of Mg and Mg + are investigated in this 
study. Prior to averaging, measurements at night or those having noise peaks in the 
investigated spectral region, caused by highly energetic particles hitting the detector 
(e.g. in the Southern Atlantic Anomaly region), are filtered from the data set. The limb 
spectra of up to 15 daily single orbits are zonally averaged before running the retrieval 
is algorithm, in order to improve the signal to noise ratio. To further reduce the noise 
in the results, all January results, February results, etc. are averaged for the entire 
measurement period. 

The algorithm retrieves densities on an altitude and latitude grid with 40 equidistant 
day side latitude intervals between 82° S and N and 80 altitude intervals with 1 km 
20 thickness between 70 and 1 50 km, additionally there are night side latitude intervals to 
overcome ambiguities of identical latitudes but different local time at high-latitudes. The 
monthly averaged results for Mg are shown in Fig. 2. 

Mg shows a peak at around 90 km with densities up to 1750 cm -3 . There is a strong 
month-to-month variability in the latitudinal distribution. The averaged yearly results are 
25 shown in Fig. 3. In the yearly average, the densities in the low latitudes are higher than 
at the high latitudes in the peak region. At the highest covered southern latitude there 
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are increased densities above 100 km, which, however, are most likely caused by the 
fact that this latitude is less covered and therefore has a larger statistical error. 

Figure 4 shows the seasonal variation of the vertical profile for low-, mid- and high- 
latitudes, and Fig. 5 shows the vertical column densities (VCDs). The month-to-month 
5 variations are large, and thus the error on the estimate of the amplitude of seasonal 
variations for mid-latitudes is relatively large. The seasonal variation of the Mg peak 
altitude is estimated to be less than 5 km. 

The VCD varies between 0.5 x 1 0 9 cm -2 and 3.5 x 1 0 9 cm -2 . The VCD between 40° N 
and 40° S is higher than for higher latitudes, which are only covered in the hemispheric 
io summer. The seasonal variation is small compared to the mean of the vertical col- 
umn densities. This relatively small variation is required to retrieve the correct mean 
value (because of non-linearities in the retrieval) as comprehensively discussed by 
Langowski et al. (2013). 


4 Seasonal variations of Mg + 

is For Mg + the same approach for data averaging is used, as described for Mg in Sect. 3. 
For Mg + there are two spectral lines, one at 279.6 nm and the second one at 280.4 nm. 
Mg + densities are independently retrieved from both lines. The differences in the peak 
region are lower than 25%. The densities on an altitude and latitude grid for the monthly 
averaged results are shown in Figs. 6 and 7. 

20 The Mg + densities peak at an altitude of 95-1 05 km with peak densities of 500- 
6000 cm -3 . The peak densities show a seasonal variation with a summer maximum 
between 25-45° in both hemispheres. In summer the highest peak altitude is roughly 
at 45° (N and S). It occurs around 105 km, which is up to 10 km higher than the lowest 
peak altitude at the equator and higher than that at the poles. There is a minimum at the 

25 equator and also a second maximum in peak altitude in the mid-latitudes on the winter 
hemisphere, that, however, is more variable than the one in the summer hemisphere. 
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The latitudinal variation of Mg + is also seen in the average over all available mea- 
surements in 2008-2012, which is shown for both lines in Fig. 8. The latitudinal depen- 
dence of the peak altitude seems to be symmetric towards the equator. However, there 
is an asymmetry in the peak density, which is close to two times as big for the Northern 
5 Hemisphere than for the Southern Hemisphere. When looking at the densities above 
the maximum density region for mid and low latitudes the density at same altitude is 
lower in the equatorial region than at 40° below 1 1 0 km, while this behavior is reversed 
above » 110 km. Furthermore, the density at the uppermost altitude level is enhanced 
in the equatorial region, pointing to the abundance of Mg + above this altitude. This may 
io be explained by enhanced vertical upward transport in the equatorial region. 

The equatorial electron distribution with transport from the equator to higher latitudes 
has been observed in ionospheric soundings. It is known as the equatorial anomaly 
or Appleton anomaly (Kendall and Windle, 1965). Typically there are fewer electrons 
at the equator than at 20° on either side of the equator due to this transport. In situ 
is satellite measurements of the metals above 120 km have been carried out using the 
Atmospheric Explorer satellite. A discussion of these data can be found in Grebowsky 
and Aikin (2002). The vertical upward transport of Mg + ions near the equator has 
been discussed, e.g. in Hanson and Sterling (1972). The sun is passing the equatorial 
region from east to west. The ionisation in the thermosphere is strong for low solar 
20 zenith angles (SZAs) which leads to denser plasmas there. The heavier ions and the 
lighter electrons gyrate in different directions in the Earth’s magnetic field pointing from 
north to south at the equator. This leads to a charge separation, which produces an 
eastwards pointed electric field and a westward current of the electrons. The magnetic 
field along the meridians and the electric field along the equator lead to an E x B-Drift 
25 perpendicular to both in the radial direction, resulting in a strong upward polarization 
field. This polarization field accelerates the ions to higher altitudes and even above the 
F-layer. However, this effect can only explain an upward transport within ±3° latitude. 
Fesen et al. (1983) showed that additionally neutral meridional winds have to be taken 
into account in order to explain vertical transport also at higher latitudes up to 30° 
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and this was experimentally shown, e.g. in Fesen and Hays (1982a, b); Gerard and 
Monfils (1 978). As for Mg the annual means show higher densities above 1 1 0 km at the 
northernmost and southernmost latitudes due to less coverage and therefore higher 
statistical errors. However, at high northern latitudes, which are cut off in the shown 
results, high densities can be retrieved when also using the stray light contaminated 
measurements. This is shown in Fig. 9. The results as well as the input raw data show 
low densities below 90 km and at different altitudes in between the high signal region, 
which makes a differentiation between a stray light effect and a true metal emission 
complicated. Should this feature be real it could be explained by the cleft ion fountain 
found by Lockwood et al., 1985, which describes the transport of charged particles 
along the magnetic field lines, which close at the pole. This transport may lift charged 
particles up to several earth radii until they become neutralized and sink down if they 
are heavy enough and are not quickly ionized again. 

Figures 10 and 11 show the seasonal variation of the vertical Mg + profile for low-, 
mid- and high-latitudes for both Mg + emission lines, and Fig. 12 shows the VCDs for 
both lines. In each of the latitude regions, a maximum in peak altitude is found for the 
summer time in the corresponding hemisphere. 

The VCD varies between 1 x 10 9 cm -2 and 11 x 10 9 cm“ 2 . A seasonal cycle with 
a summer maximum is observed and the strongest variations are observed between 
20° and 40°. Figure 13 shows the ratio of Mg + and Mg VCDs. For low- and mid- latitudes 
as well as the high latitudes in the southern hemispheric summer the ratio is in between 
0.5 and 4. The seasonal variability of the ratio follows the one of Mg + , as Mg does not 
show a strong seasonal variability. 


5 Comparison to other measurements 

Mg and Mg + cannot be observed from the ground, because the wavelengths of the lines 
are below 300 nm and the emission signal is absorbed strongly by ozone in the strato- 
sphere. As a result, only few measurements are available and rely on observations from 
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sounding rockets and satellites. The first rocket-borne ion mass spectrometer measure- 
ment of metal ions was performed by Johnson and Meadows (1955) (May 1954, White 
Sands) and enhanced ion signals between 93 and 124 km were found. According to 
Grebowsky and Aikin (2002) approximately 50 flights of rocket-borne mass spectrome- 
5 ters, probing the region between 80 and 130 km have been performed since then (until 
2002). Other density profiles for Mg + can be found in Istomin (1963); Narcisi and Bai- 
ley (1965); Narcisi (1971); Aikin and Goldberg (1973); Philbrick et al. (1973); Zbinden 
et al. (1975); Herrmann et al. (1978); Kopp and Herrmann (1984); Kopp et al. (1985b, 
a); Kopp (1997); Roddy et al. (2004). Date, local time, latitude and reference publica- 
io tions of these and other rocket flights can also be found in Table 1 of Grebowsky et al. 
(1998). 

A disadvantage for the comparison with the new dataset presented is that most of 
these measurements were performed during special events with E s layers, aurora, me- 
teor showers, stratospheric warmings or noctilucent clouds present. Often sporadic 
is Mg + layers occured between 105-1 10 km and/or 120 km. The Mg + peak altitude in 
most of these measurements can be found between 90 and 95 km altitude and the full 
width at half maximum (FWHM) of the layers is on the order of 5-1 0 km, but sometimes 
the FWHM is only 1 km. 

When comparing in-situ mass spectrometer measurements with satellite remote 
20 sensing results, it has to be remembered that the in-situ measurements are much more 
localized and limited to the direct vicinity of the rocket. In contrast, remote sensing tech- 
niques typically cover a large volume with horizontal distances along and perpendicular 
to the viewing direction of several hundreds of kilometers. This results in smoother lay- 
ers with a larger FWHM for the remote sensing method. The width of the Mg + layers 
25 as well as the peak densities presented in this study are in good agreement with the 
in-situ rocket measurements. However, the strong latitudinal dependence of the peak 
altitude observed in Figs. 10 and 11 is not found in the in-situ rocket data. In addi- 
tion to the in-situ measurements with rockets, there are also airglow measurements 
available, from rockets, space shuttles and satellites. However, to retrieve density in- 
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formation from this method radiative transfer models as well as inversion techniques 
and computational power is needed. Therefore, often only slant column information is 
available (see, e.g. Boksenberg and Gerard, 1973 or Gerard and Monfils, 1974). In An- 
derson and Barth (1971) (Summer, « 40° N) the region up to 1 06 km altitude is scanned 
5 during a sporadic E s layer event, and the peak altitude of the Mg + is not reached dur- 
ing this flight, i.e. was higher than 106 km. No Mg signal above the instrumental noise 
was observed in these spectra. The region above the peak altitude from 150 km up to 
the F-layer and above is investigated in Gerard and Monfils (1978); Fesen and Hays 
(1 982a); Mende et al. (1 985); Gardner et al. (1 995, 1 999) and typically shows less than 
io 100cm -3 Mg + ions at 150 km. This is in good agreement with the profiles described 
in Sect. 4. In some plots higher densities can be seen in the profiles retrieved from 
the limb observations of SCIAMACHY. These are explained by retrieval artifacts on the 
edge of the retrieval grid and vertical constraints tuned for the main peak. These allow 
small oscillations, compared to the main peak, in regions with lower densities. Since 
is these artifacts appear at the upper edge of the profile it also hints that there is a signif- 
icant density above the highest tangent altitude for Mg + . In Minschwaner et al. (2007) 
a combined NO and Mg + retrieval for satellite limb measurements is shown, which is 
in good agreement with the results in Sect. 4 in terms of the Mg + concentrations at 
peak altitude and at the upper edge of the profile. However, the Mg + peak altitude is at 
20 90 km. Taking into account the coarser sampling (every 7 km), higher statistical errors 
and a tangent height offset (±4 km) in Minschwaner et al. (2007), the agreement is 
reasonable. Tangent height offset of a similar magnitude was also an issue for SCIA- 
MACHY (von Savigny et al., 2005), however, this error was reduced to (±200 m) (von 
Savigny et al., 2009). The NO Band emission, which overlaps with the Mg + lines and 
25 the Mg line, is on the same order of magnitude as the Mg + lines in Minschwaner et al. 
(2007) and even bigger than the Mg emission, which made a NO correction neces- 
sary. The SCIAMACHY MLT dataset does not show these strong NO lines at 280 and 
285 nm. NO in this region is very sensitive to solar activity. Results for NO retrievals 
from the same SCIAMACHY dataset used in this study are reported by Bender et al. 
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(2013). These exhibit densities of the same order of magnitude as the equatorial den- 
sity plot for NO in Minschwaner et al. (2007), only at high-latitudes in summer and 
winter and from late 201 1 to 2012. However, we did not observe clear NO signals in 
the vicinity of the Mg/Mg + lines in this period. A time series of Mg + vertical columns 
s covering several years and retrieved from SBUV nadir measurements was presented 
by Joiner and Aikin (1996). These measurements were performed approximately 1 day 
per month, with a spectral resolution of 1.13 nm and a spectral sampling every 0.2 nm 
(a factor 5 for the resolution, and 2 for the sampling, worse than SCIAMACHY). The 
results in Joiner and Aikin (1996) are in very good agreement with the results obtained 
io in this study, especially when comparing Figs. 1 0 and 1 1 in Joiner and Aikin (1 996) with 
Fig. 12 in this study. Mg was also investigated in Joiner and Aikin (1996). However, the 
average VCD for the rare profiles, where the signal was significant, is 4 x 10 10 cm -2 , 
which is a factor 10 more than the VCDs in Fig. 5. 

Satellite measurements with long time series and daily coverage are available from 
is GOME and GOME-2 in nadir mode and from SCIAMACHY in nadir and nominal limb 
mode. All 3 instruments have a similar spectral resolution. The GOME dataset was 
investigated for Mg + and Mg VCDs in Correira et al. (2008, 2010) and Correira (2009). 
Figures 1 and 2 in Correira et al. (2008) show Mg and Mg + VCDs for 1996 and 1997 as 
well as the Mg + to Mg ratio for latitude intervals from 0-10° and 30-40° for both hemi- 
20 spheres. For the low-latitudes, where there is less seasonality the VCDs for Mg + are 
at 6-7x10 9 cnrf 2 which is higher than in Fig. 12, the Mg densities are at 3x10 9 cnrT 2 
which is a similar magnitude to that in Fig. 5. In the equatorial region the Mg columns 
agree better, and a higher VCD in the Nadir results for Mg + can be explained by the 
thermospheric part of Mg + which is not part of the VCDs in Fig. 12. At mid-latitudes 
25 Mg + shows a strong seasonal cycle with a summer maximum in Correira et al. (2008). 
This seasonality is quite symmetric for both hemispheres and the bigger VCDs can 
be found in the Southern Hemisphere. Furthermore, the summer maximum at mid lat- 
itudes does not exceed the VCDs at low-latitudes. For Mg a summer maximum is ob- 
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served which is more pronounced in the Southern Hemisphere. The Mg + to Mg VCD 
ratios in Correira et al. (2008) are in good agreement with our results shown in Fig. 13. 

The nominal limb mode dataset, covering the tangent height range from the surface 
up to 92 or 105 km (only until early 2003), respectively, combined with the nadir mode 
5 dataset from SCIAMACHY from 2002 to 2007 was investigated in Scharringhausen 
et al. (2008b, a) and Scharringhausen (2007). The peak region for the ions could not 
be fully resolved with the nominal limb mode, however, the agreement for ions is quite 
good showing only slightly higher VCD in Scharringhausen (2007) compared to the 
results in this study. In the Southern Hemisphere the seasonal cycle was not identi- 
io fied well in Scharringhausen (2007), which agrees with the results presented in Fig. 12 
showing a weaker seasonal cycle in this region. At the very beginning of the Scharring- 
hausen VCD dataset in 2002 (e.g. Figs. 69 and 70 in Scharringhausen, 2007) slightly 
higher VCDs for both species are observed than for the other years. This may come 
from increased NO signals near the Mg/Mg + lines during solar maximum. Another 
is reason for the higher Mg + VCDs in 2002 may also be the change of the maximum 
tangent altitude from 105 to 95 km. More differences to Scharringhausen (2007) can 
be found in the Mg data, since the Mg line is much more affected by newly added ra- 
diative transfer features. There was no correction of the Ring effect in Scharringhausen 
(2007), which led to high densities below 80 km and a density maximum at the lower 
20 edge of the retrieval boundaries at 70 km. This inelastic scattering contribution addi- 
tionally adds a seasonal variation to the dataset. Furthermore, the Mg line at 285.2 nm 
is much more affected by self absorption of the emission, which was not considered in 
Scharringhausen (2007), so the Mg VCDs were smaller than reported here. However, 
most of the findings from the further analysis of the Scharringhausen dataset are only 
25 weakly affected by these differences and are still valid. 

In summary it can be concluded, that the Mg/Mg + results presented here are in good 
agreement with previous satellite and rocket instruments. The most striking difference 
to previous measurements is the strong latitudinal dependence of the peak altitude of 
Mg + with differences of up to 10 km for different latitudinal zones. 
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6 Comparison of Mg and Mg + observed by SCIAMACHY and modelled by the 
WACCM model 

6.1 WACCM Model 

WACCM is a high-top chemical-dynamical model, which simulates the altitudes from 
the earth’s surface up to 140 km (Garcia et al., 2007; Marsh et al., 2007, 2013b). Here 
we use the specified dynamics SD-WACCM, which is nudged by the GEOS5 meteoro- 
logical dataset (including temperature, specific humidity, horizontal winds) below 60 km. 
This is the same version used by Feng et al. (2013) and Marsh et al. (2013a) who 
successfully incorporated iron and sodium chemistry into WACCM. Here we develop 
a global model of meteoric magnesium in the atmosphere by combining three compo- 
nents: a treatment of the injection of meteoric constituents into the atmosphere (Feng 
et al., 2013; Marsh et al., 2013a), a description of the neutral and ion-molecule chem- 
istry of magnesium in the mesosphere and lower thermosphere (Whalley and Plane, 
2010; Plane and Whalley, 2012) and WACCM. The Mg meteoroid injection function 
(MIF) has a similar seasonal variability as those of Fe and Na. The minimum MIF is 
in the spring with an ablation flux of 1300 atoms cm _2 s _1 , and the maximum occurs 
in the autumn with 2400cm _2 s _1 . The peak ablation height is around 95 km. These 
fluxes are obtained by scaling the Mg MIF to match the observed Mg column density. 
This requires dividing the Mg ablation flux by a factor of 15, relative to that of Na (Marsh 
et al., 2013a). Since the relative chondritic abundance of Mg is 16 times that of Na, this 
means that the ablation fluxes of Mg and Na are similar in the model. Tables 1-3 list 
the Mg chemical reactions and their rate coefficients added to the standard chemistry 
scheme in this study. 

An important magnesium reservoir on the underside of the Mg layer is Mg(OH) 2 . 
This is reduced back to Mg via MgOH by reaction with H atoms (R12). The subsequent 
reaction with H atoms (R12b) is most likely faster than R12, and so the reaction rate 
is not explicitly in Table 2 (Plane and Whalley, 2012). The polymerization of Mg0 2 H 2 
to form meteoric smoke is parameterised by a dimerization rate coefficient (reaction 

16 


if) 

O 

c 

if) 

c n 

o' 

-o 

"0 

0) 

"O 

CD 


if) 

o 

c 

if) 

if) 

o' 

o 

"□ 

0) 

■O 

CD 


ACPD 

14, 1-49, 2014 


Global investigation 
of the Mg atom and 
ion layers 

M. Langowski et al. 


Title Page 


Abstract ■ Introduction 


Conclusions ■ References 


if) 

o 

c 

if) 

c n 

o' 

o 

~0 

0) 

"O 

CD 


Tables 


Figures 



D 

in' 

o 

c 

in 

in 

o' 


"0 


0) 

“O 

CD 


Full Screen / Esc 

Printer-friendly Version 
Interactive Discussion 


















R13 in Table 2). This reaction should be essentially at the high pressure limit even in 
the upper mesosphere because of the large number of atoms in the dimer and the 
large binding energy (268 kJ mol -1 , calculated at the B3 LY P/6-31 1 +g(2d, p) level of 
electronic structure theory (see, e.g. Frisch et al., 2009). The capture rate coefficient is 
then increased to 9 x 1 0 -1 °cm 3 molecule -1 s -1 to account for the concentration of other 
metallic species (e.g. NaHC0 3 , FeOH etc) with which Mg0 2 H 2 can also polymerize. 

The model was run for the period 2004 to 201 1 , when the GEOS5 analysis data are 
available. To derive the modeled climatologies of temperature, Mg and other chemical 
constituents, we use the model output from 2005 to 201 1 . 

6.2 Comparison of WACCM model and SCIAMACHY measurements results 

Figure 14 shows the monthly mean results for 7yr of WACCM simulations of Mg from 
2005 to 201 1 and Fig. 15 shows the same for Mg + . Figure 16 shows the VCDs for both 
Mg and Mg + . Mg shows a clear seasonal cycle with a winter maximum which is most 
pronounced at the poles. The peak altitude is nearly constant, but small seasonal vari- 
ations can be found. E.g. in February the peak altitude at the poles is roughly 85 km, 
while it is [5] km higher at 90 km in the equatorial region. Mg + shows a seasonal cycle 
with a summer maximum. In addition, there is an increased Mg + density at the poles, 
even in the winter hemisphere. The Mg + peak altitude is close to 95 km and shows no 
strong variation, with latitude and time. The seasonal variations in the VCD profiles for 
Mg and Mg + are very similar to Na and Na + profiles (Fig. 3 in Marsh et al., 2013a). 
In contrast to the SCIAMACHY measurements, which are made at one particular lo- 
cal time, daily averaged output is used for WACCM. We also colocated the WACCM 
dataset to SCIAMACHY local time and latitudinal coverage, but found only small differ- 
ences to the non collocated daily average mean of WACCM. There is a diurnal variation 
of the vertical column density in the WACCM data, but nearly no diurnal variation for 
the the vertical profile shape. Figures 17 and 18 show the direct comparison of SCIA- 
MACHY and WACCM annual means for Mg and Mg + . We used the colocated data for 
both plots, since these are the only plots with significant differences between uncollo- 
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cated daily means and colocated data, for latitudes higher than 40°. Figures 19 and 20 
show the comparison for the vertical profile only. 

In the annual mean profile the peak density of Mg for WACOM is higher than for 
SCIAMACHY. However, the peak width of the SCIAMACHY profile is wider than for 
5 WACOM and the VCDs are very similar, which is expected, because the Mg MIF is 
scaled (see above). The peak altitude of WACCM (88 km) and SCIAMACHY (90 km) 
results agree very well considering the 3.3 km step size in the SCIAMACHY measure- 
ments and the vertical sampling step size of 2 km for WACCM and 1 km for SCIA- 
MACHY. One explanation for the wider SCIAMACHY peak is a strong Ring effect in- 
io fluence (Grainger and Ring, 1962) for the Mg spectral line, which has not adequately 
been corrected. At the lower and upper peak edge, the SCIAMACHY number den- 
sities are not zero. This may be related to the requirement in the retrievals that the 
densities must be positive, which is needed for stability in the retrieval. However, the 
densities are smaller than the single day error there. The agreement of WACCM and 
is SCIAMACHY is better for Mg + . The WACCM Mg + peak and column densities are of 
similar magnitude as the SCIAMACHY results and also show the same seasonal vari- 
ability with a summer maximum. There are also small differences: the maximum in the 
WACCM data is not reduced when passing the equatorial region during equinox and 
also has a stronger extension to the higher latitudes in summer. Furthermore, a sec- 
20 ond maximum appears in the winter hemisphere at high-latitudes, which is not covered 
by SCIAMACHY measurements and so is not observed. The WACCM results also do 
not show the strong seasonality and latitude dependence of the peak altitude, that can 
be found in the SCIAMACHY results. The peak altitude at 95 km is in slightly better 
agreement with the rocket data than with SCIAMACHY. One reason for some of these 
25 discrepancies between model and observations is that WACCM does not include the 
Lorentz v x£-force, as discussed in Feng et al. (201 3). In a global average of the profile, 
however, the differences are small. The lower peak edge of Mg + is in better agreement 
than for Mg, since, e.g. the Ring effect influence is much weaker for Mg + than for Mg. 
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The agreement of the peak edge becomes worse with larger distance from the peak 
altitude. Nevertheless, the FWHM of the peak is in good agreement. 


7 Summary and conclusions 

Monthly averaged global vertical and latitudinal density distributions of Mg and Mg + 
retrieved from the SCIAMACHY limb MLT measurements from 2008 to 2012 have been 
shown. The MLT measurements cover the vertical region between 70 km and 150 km 
and have a good resolution in the peak region. Mg shows highly variable results and 
variabilities are mostly affected by the measurements error and no clear seasonal cycle 
can be observed for Mg. The peak altitude of Mg is nearly constant at » 90 km for all 
latitudes and times. Mg + shows a clear seasonal cycle with a summer maximum in the 
peak density, most pronounced at mid-latitudes. Furthermore, Mg + shows a seasonal 
variation of the peak altitude, with higher altitudes in the summer, as well as a latitudinal 
dependence of the peak altitudes with up to 5-1 0 km higher altitudes at mid-latitudes 
compared to the equatorial peak altitude at 95 km. 

The SCIAMACHY results are reasonably consistent with previously published rocket 
profiles. There is a good agreement for the peak density values as well as the seasonal- 
ity for Mg + . However, the rocket measurements do not show the latitudinal dependence 
of the peak altitude, but typically show peak altitudes at slightly below 95 km. 

The SCIAMACHY results were compared with WACCM model results. There is a dis- 
agreement in the seasonal variation for Mg between SCIAMACHY and WACCM results, 
with WACCM showing a clear seasonal cycle with a winter maximum at the poles. Fur- 
thermore, the WACCM peak densities are roughly a factor 2 bigger than the SCIA- 
MACHY peak densities. The peak altitude and peak shape are in good agreement with 
SCIAMACHY, showing a slightly wider profile. The combination of higher peak densi- 
ties for WACCM, but wider peak profile for SCIAMACHY leads to very similar values 
for the VCDs. The agreement between WACCM and SCIAMACHY is better for Mg + 
than for Mg. Both SCIAMACHY and WACCM show a clear seasonal cycle with a sum- 

19 


if) 

O 

c 

if) 

c n 

o' 

-o 

"0 

0) 

"O 

CD 


if) 

o 

c 

if) 

if) 

o' 

o 

"□ 

0) 

■O 

CD 


ACPD 

14, 1-49, 2014 


Global investigation 
of the Mg atom and 
ion layers 

M. Langowski et al. 


Title Page 


Abstract ■ Introduction 


Conclusions ■ References 


if) 

o 

c 

if) 

c n 

o' 

o 

~0 

0) 

"O 

CD 


Tables 


Figures 



D 

in' 

o 

c 

in 

in 

o' 


"0 


0) 

“O 

CD 


Full Screen / Esc 

Printer-friendly Version 
Interactive Discussion 


















mer maximum and the peak densities are in good agreement. However, like the rocket 
measurements, WACOM does not show a latitudinal dependence of the peak altitude. 

The peak shape of Mg + agrees very well for WACCM and SCIAMACHY in the high 
density region, where the SCIAMACHY error is low. 

a? 
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Table 1. Magnesium Photolysis Chemistry added into WACCM. 


No. 

Reaction 

Rate/s 1 

Reference 

J1 

Mg + hv — * Mg + + e _ 

y'i =4x1 0 -7 

Swider (1984) 
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Table 2. Neutral chemistry of magnesium added into WACOM. 


No. 

Reaction 

Rate/cm 3 molecule 1 s 1 

Reference 

R1 

Mg + 0 3 — » MgO + 0 2 

k : =2.3x 10“ 10 exp(-139 /T) 

Plane and Helmer (1995) 

R2 

MgO + O — » Mg + 0 2 

k 2 = 6.2 x 10' 10 x (T/295) 0 - 167 

Whalley and Plane (2010) 

R3 

MgO + 0 3 — ► Mg0 2 + 0 2 

k 3 = 2.2 x 10 _1 °exp(-548/T) 

Plane and Helmer (1995) 

R4 

Mg0 2 + O — » MgO + 0 2 

k 4 = 7.9 x 10" 11 exp(7/295) 0,167 

Whalley and Plane (2010) 

R5 

MgO + H 2 0 + M -» Mg0 2 H 2 + M 

k 5 = 1.1 x 10- 26 x (7/200)" 1 - 59 

Rollason and Plane (2001) 

R6 

Mg0 3 + H 2 0 — » Mg0 2 H 2 + 0 2 

k 6 = 1.0 x 10“ 12 

Rollason and Plane (2001) 

R7 

MgO + 0 2 + M — > Mg0 3 + M 

k 7 = 3.8 x 10 -29 x (7/200)" 1 ' 59 

Rollason and Plane (2001) 

R8 

MgO + C0 2 + M — > MgC0 3 + M 

k 8 = 5.9 x 10 -29 x (7/200) -0 ' 86 

Rollason and Plane (2001) 

R9 

MgC0 3 + O — » Mg0 2 + C0 2 

k 9 = 6.7x 10“ 12 

Plane and Whalley (2012) 

R10 

Mg0 2 + 0 2 + M — > Mg0 4 + M 

/r 10 = 1 .8 x 10 -26 x (7/200) -2 ' 5 

Plane and Whalley (2012) 

R11 

Mg0 4 + O -» Mg0 3 + 0 2 

k u = 8.0 x 10~ 14 

Plane and Whalley (2012) 

R12 

Mg0 2 H 2 + H — » MgOH + H 2 0 

/r 12 = 1 .0 x 1 0" 1 1 x exp(-600/7) 

Plane and Whalley (2012) 

R12b 

MgOH + H — » Mg + H 2 0 

faster than R12 

see text 

R13 

Mg0 2 H 2 + Mg0 2 H 2 — » (Mg0 2 H 2 ) 2 

/r 13 = 9.0 x 10~ 10 

see text 

R14 

Mg0 3 + H — > MgOH + 0 2 

/r 14 = 2.0 x 10~ 12 

Plane and Whalley (2012) 
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Table 3. Ion-molecule chemistry of magnesium added into WACCM. 


No. 

Reaction 

Rate/cm 3 molecule 1 s 1 

Reference 

R15 

Mg + + 0 3 — » MgO + + 0 2 

/r 15 = 1 .2 x 10 -9 

Whalley and Plane (2010) 

R16 

MgO + + O — » Mg + + 0 2 

/r 16 = 5.9 x 10 _1 ° 

Whalley and Plane (2010) 

R17 

MgO + O 3 — > Mg + 2 O 2 

/r 17 = 1 .8 x 10 _1 ° 

Whalley and Plane (2010) 

R18 

MgO + + 0 3 — ► Mg0 2 + 0 2 

/r 18 = 3.3 x 10 _1 ° 

Whalley and Plane (2010) 

R19 

Mg + + N 2 + M — » MgN 2 + + M 

/r 19 = 1 .8 x 1 0 _ 30 (r / 200) -1 ' 72 

Whalley and Plane (2010) 

R20 

Mg + + 0 2 + M — > MgN 2 + + M 

k 20 = 2.4 x 1(T 30 (7"/200 ) -1 ' 65 

Whalley and Plane (2010) 

R21 

MgN^ + e — » Mg + N 2 

k 2: = 3.0 x icr 7 (r/ 2 oo) -0 ' 5 

Plane and Whalley (2012) 

R22 

MgCO^ + e — > Mg + C0 2 

k 22 = 3.0 x 10 -7 (r/200) -0 ' 5 

Plane and Whalley (2012) 

R23 

MgH 2 0 + + e — » Mg + H 2 0 

k 23 = 3.0 x 10 -7 (r/200) -0 ' 5 

Plane and Whalley (2012) 

R24 

MgO + + e — > Mg + O 

k 24 = 3.0 x icr 7 (r/ 2 oo) -0 ' 5 

Plane and Whalley (2012) 

R25 

Mg + + e“ — > Mg + hv 

k 25 = I.OxlO -12 

Plane and Helmer (1995) 

R26 

Mg + 0 2 + — » Mg + + 0 2 

k 2e = 1 .2 x 1 cr 9 

Anderson and Barth (1971) 

R27 

Mg + N 2 + — > Mg + + N 2 

k 27 = 8.2 x 1CT 10 

Anderson and Barth (1971) 

R28 

MgN 2 + + 0 2 — » Mg0 2 + + N 2 

k 28 = 3.5 x 10 _12 

Rollason and Plane (2001) 

R29 

Mg0 2 + + O — > MgO + + 0 2 

k 29 = 6.5 x 1CT 10 

Rollason and Plane (2001) 

R30 

MgO + O 3 — ► Mg + 20 2 

k 30 = 1 .8 x 10 _1 ° 

Whalley and Plane (2010) 

R31 

Mg + + H 2 0 + M -> MgH 2 0 + + M 

/r 31 = 2.3 x 1 0 _28 (7" /200) -2 ' 53 

Martinez-Nunez et al. (2010) 

R32 

Mg + + C0 2 + M — » MgC0 2 + + M 

k 32 = 4.6 x 10 -29 (7"/200) -1 ' 42 

Whalley and Plane (2010) 
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Fig. 1. Local time of SCIAMACHY measurements at different latitudes for different seasons of 
the year. 
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Fig. 2. Latitudinal variation of monthly (averaged over all January, February, etc. results from 
2008-2012) and zonally averaged Mg density profiles retrieved from the 285.2 nm line. The 
Mg layer peaks around 90 km altitude and has a FWHM of about 15 km. Mg shows strong 
month-to-month variations, but there is no clear seasonal cycle to be observed. 
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Fig. 3. Latitudinal variation of Mg profiles averaged over all available data for the years 2008- 
2012. Note that only between 45° N and S measurements from all 12 months contribute to 
the averages, while for high-latitudes there is only summer coverage and/or near-terminator 
measurements. 
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Fig. 4. Seasonal variation of the vertical Mg profile for different latitudinal zones (left: Northern 
Hemisphere, right: Southern Hemisphere, top: low-latitudes, mid: mid-latitudes, bottom: high- 
latitudes). The peak altitude is at 90 km for all latitudes with variations of ±5 km. The month-to- 
month variations in peak altitude and density are bigger than any seasonal variations. 
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Fig. 5. Seasonal variation of the Mg vertical column density (VCD) between 70 km and 1 50 km 
for different latitudinal regions. The VCD varies between 0.5 x 1 0 9 cm -2 and 3.5 x 1 0 9 cm -2 . No 
clear seasonal cycle is observed. Variations between consecutive days with measurements (at 
least 14 days difference) are even higher than the inter-monthly variations. 
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Fig. 6. Latitudinal variation of monthly and zonally averaged Mg + profiles (for the years 
2008-2012) retrieved from the Mg + line at 279.6 nm. Mg + shows a seasonal cycle with a sum- 
mer maximum, which is especially pronounced in the region between 25° and 45° latitude in 
the summer hemisphere. Furthermore, in this region the peak altitude is about 105 km, which 
is up to 1 0 km higher than at the equator or at the poles. 
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Fig. 7. Latitudinal variation of monthly and zonally averaged Mg + profiles (for the years 
2008-2012) retrieved from the Mg + line at 280.4 nm. The results from the 280.4 nm line are 
in good agreement with the results based on the 279.6 nm line shown in Fig. 6. The 280.4 nm 
line retrievals show slightly higher densities above 110 km and a slightly smaller peak den- 
sity, especially when the peak density is high. Note that the peak at 150 km and high northern 
latitude is probably an artifact of the retrieval. 
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Fig. 8. Latitudinal variation of Mg + profiles averaged over all available data for 2008-2012 for 
the Mg + line at 279.6 nm (left) and the Mg + line at 280.4 nm (right). Results for both lines show 
an asymmetric behaviour between the Northern Hemisphere and the Southern Hemisphere 
with a more pronounced summer maximum in the Northern Hemisphere. 
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Fig. 9. Mg + profiles averaged over all available data for 2008-201 2 for the Mg + line at 279.6 nm 
(left) and the Mg + line at 280.4 nm (right). The stray light contaminated measurements at the 
start of each orbit are also used here extending the coverage to higher northern latitudes and 
even the ascending node side. High Mg + densities are observed at the north pole. 
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Fig. 10. Seasonal variation of the vertical profile of Mg + (279.6 nm line) for different latitudinal 
zones (left: Northern Hemisphere, right: Southern Hemisphere, top: low-latitudes, mid: mid- 
latitudes, bottom: high-latitudes), averaged over all available observations in 2008-2012. For 
all these latitudinal zones a seasonal cycle in the peak altitude with a summer maximum and 
a winter minimum can be observed. The seasonal variations are similar or larger than the 
month-to-month variations. The maximum peak altitude is 5-1 0 km higher than the minimum 
peak altitude. For mid-latitudes, the peak altitude in the winter hemisphere is still higher than 
the peak altitude at the equator (see also Fig. 6). 
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Month 


Fig. 11. Seasonal variation of the vertical profile of Mg + (280.4 nm line) for different latitudinal 
zones (left: Northern Hemisphere, right: Southern Hemisphere, top: low-latitudes, mid: mid- 
latitudes, bottom: high-latitudes), averaged over all available observations in 2008-2012. The 
results are in good agreement with the results in Fig. 10. 
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Fig. 12. Seasonal variation of the Mg + vertical column density (VCD) between 70 km and 
150 km and for different latitudinal regions, retrieved from the 279.6 nm line (left) and the 
280.4 nm line (right). A clear seasonal cycle with a summer maximum can be observed. The 
summer maximum in the Northern Hemisphere has higher values compared to the one in the 
Southern Hemisphere. The highest variability can be found between 25° and 50° in both hemi- 
spheres. 
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Fig. 13. Seasonal variation of the Mg + (left 279.6 nm, right 280.4 nm) to Mg ratio for different 
latitudes. At high northern latitudes the ratio may exceed the value of 6. For mid-latitudes the 
ratio varies between 0.5 and 4 and shows a summer maximum, which is in good agreement 
with Correira et al. (2008). 
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Fig. 14. Monthly mean WACOM results for Mg. The model results show a clear seasonal cycle 
with a winter maximum, most pronounced in the high-latitudes. 
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Fig. 15. Monthly mean WACOM results for Mg + . The model results show a clear seasonal 
cycle with a summer maximum, most pronounced in the high-latitudes. There is also a second 
smaller maximum at high-latitudes for the winter hemisphere. 
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Fig. 16. VCD timelines for different latitudes for Mg (left) and Mg + (right). Mg VCDs vary be- 
tween 0.5 x 1 0 9 cm -2 and 8 x 1 0 9 cm -2 and show a clear seasonal cycle with a winter maximum, 
which is most pronounced at the poles. Mg + VCDs vary from 1 xIO 9 cm -2 to 5 x 10 9 cm -2 , with 
most VCDs between 2 x 10 9 cnrf 2 and 3 x 10 9 cm -2 . A clear seasonal cycle with a summer 
maximum for Mg + is observed. SCIAMACHY results of VCDs are shown in Figs. 5 and 12. 
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Fig. 17. Annual means of Mg for SCIAMACHY (left) and WACCM (right). SCIAMACHY shows 
(roughly a factor 2) smaller densities and a wider vertical profile. The averaging for WACCM 
considers only the time periods, where SCIAMACHY measurements are available. 
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Fig. 18. Annual means of Mg + for SCIAMACHY (left, 279.6 nm) and WACCM (right). SCIA- 
MACHY shows a wider vertical profile, and a stronger latitudinal dependence of the peak al- 
titude. The peak densities are in very good agreement. The averaging for WACCM considers 
only the time periods, where SCIAMACHY measurements are available. 
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Fig. 19. Global annual mean vertical profile of Mg for SCIAMACHY and WACOM. The right 
plot shows the profiles normalized to the maximum of the WACOM profile, to better compare 
the shape of the profiles. The smoothed profile is smoothed with a triangular function with 
a base width of 4 km. The VCD between 70 and 140 km is 2.52 x 10 9 cm“ 2 for SCIAMACHY 
and 2.49 x 1 0 9 cm" 2 for WACOM. 
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Fig. 20. Vertical profiles of Mg + for SCIAMACHY and WACOM at the equator. The right plot 
shows the profiles normalized to the maximum of the WACCM profile, to better compare the 
shape of the profiles. The profile shape for the 279.6 nm line matches the WACCM profile shape 
very well. The VCD between 70 and 1 40 km at the equator is 3. 1 x 1 0 9 cm -2 and 2.8 x 1 0 9 cm -2 
for SCIAMACHY (279.6/280.4 nm) and 2.4 x 10 9 cm -2 for WACCM. 
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